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f***) Lattice effects on the kink families of two models for one-dimensional nonlinear Klein-Gordon sys- 

tems with double-well on-site potentials are considered. The analytical expression of the generalized 
Peierls-Nabarro pinning potential is obtained and confronted with numerical simulations. 
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Many physical problems have greatly contributed to the interest on double-well potential© WP) models Q. Among 
, the most important ones are the mechanism of protonic conductivity in molecular systems [^],p|, structural instabilities 
in one-dimensional(lD) Perovskyte-like systems Q as well as Hydrogen-bonded ferroelectrics || , and Peierls transition 
in ID conductors [pi. The first problem ows its interest to the fundamental role of proton migrations among the 
surrounding (heavy) ions of the ID molecular system. The ID character of these systems sits in the very high 
conductivity along the Hydrogen-Bond(HB)-about 10 3 times greater-compared to the two other directions Q. This 
strong anisotropic conductivity relies on structural defects, that is ionic and bonding defects Q . Defects are consequent 
upon proton hopping and bond distortions across the energy barrier between the two most stretched proton positions 
along the HB. When evoking ferroelectricity and the Peierls transition(an instability of the electronic structure of 
a ID electron-phonon system), the most natural framework is that of the Landau theory of second-order phase 
[ transition The "DWP" scheme in this theory is appropriate for such phenomena as the proton ordering transitions 
'55 ' ^ n ^ e P seuc ' X) ~lI- > ferroelectrics KH2PO4 || and structural instabilities in ID organic and inorganic conductors [[|. 
For all these systems, in an idealized context in the absence of external bias(electric or magnetic fields), defect 
nucleations(soliton lattice) |n| and the intrinsic discrete structure of the background substrate lattice |ll[] are the 
O |' two main processes influencing conductivity. In general they give rise to an activated transport property that consists 
either in a conducting phase or in an insulating phase |L^|L3| ]. In concrete terms, by the Kramers theory of diffusion it 
' is easy to point out that the kink diffusion coefficient in a discrete medium will be weighted by an Arhenius-like factor 
R(lu p ) which accounts for the escape rate of kink from the lattice trapping. This factor depends on the vibration 
\Q • frequency u p of the kink in its attempts to jump across the pinning potential barrier, thereby determining its lifetime 
ON I in the discrete medium. At low temperatures, the discrete kink diffusion coefficient is of the form 
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D ~ D R(u p ), R(u„) ~ u p eM~PE PN ) (1) 

Epn in this formula is the lattice-induced trapping potential barrier, the kink pinning frequency lj p is directly con- 
nected to this physical quantity. The activated diffusion coefficient given above is not only interesting in its own right, 
& , but also it furnishes relevant knowledge for capturing the essential features of all other transport parameters namely, 
the kink mobility and conductivity. 



While DWP models have so become useful in the understanding of physical problems, insights gained from the earliest 



model i.e. the 4 rapidly saturated owing to several shortcomings. For instance, the rigidity prevents from getting 
further in the theory as to analytically account for the possibility of shifting the potential minima and hump to adjust 
the model to a desired context. This has been particularly dramatic for molecular solids which backbones are often 
so flexible that chemical bonds must continuously relax or contract to avoid bond breakings. To overcome these 
shortcomings, the Morse and double-Morse(DM) potentials were considered. Thus, the DM potential has been used 
in the recent studies of the protonic conductivity Jl4[] and the proton ordering transitions p5[ . For this last context, 
it has even been argued |l5| that the DM model is to date the best candidate giving relatively good account of the 
so-called "geometric effect" . 



In addition to the DM, several other parametrized DWP models are currently present in the literature [|16 17 1. By 
their essential virtues they allow theoretical manipulations at one wish. The present work deals with two versions 
among the most general ones. A main question in the theory of parametrized DWP (PDWP) is their analytical 
tractability. Indeed, except the soliton solutions that are often easy to obtain, several other physical quantities are 



explicitely unaccessible. Even their statistical mechanics are most often proceeded numerically 18 , or by transfer 
integral methods which however can open only very narrow paths in the rather wide range of thermodynamic pa- 
rameters of the system For the present two variants we will discuss the lattice discreteness properties to both 
analytical and numerical viewpoints. We start by deriving a generalized analytical expression of the Peierls-Nabarro 
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potential(PNP) [l2|| . Next, we carry out a numerical analysis to get into the intimate discrete structure of each model. 
By this last way we also provide an adequate framework to check the exactness of our analytical results. 
Starting, consider a ID nonlinear Klein-Gordon system whose dynamics can be described by the Hamiltonian: 



H = M Y J 



n=l 
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V(<j>, fi) is the PDWP assumed of the form [p 
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W) = o I — sinh"(a0)-l 



The quantities a and a can take different values, namely: 



a = a 



a = fi 



For this value set, (|^) is a DWP with variable positions of degenerate minima but fixed potential hump. 
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(1 + fi 2 )arsinh 2 (fj,)' 



a = arsinh(fi) 
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for this second value set (|^) is a DWP with fixed minima but variable potential hump. In the continuum limit(x = n£, 
where I is the lattice spacing), the PDWP admits single-kink solutions explicitly given by: 



Where we set 7 



4> n (s = nt — vt, u) = ±—artanh 



and call d(p) the kink width: 
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dk in this last relation is the kink "bare" width. At this step it is instructive to remark that according to the two 
model parameters, only the first will give rise to kink family whose widths effectively vary(decrease) as function of 
(with increasing) /i, whereas all kinks of the second family have their widths bounded to the constant dk whatever //. 
However, their "kink shape" is deeply affected by the deformability parameter as noticed elsewhere fll7|| . 
As our primary goal is the kink-lattice interactions, proceeding with we appeal to the PNP approach which assumes 
the existence of an effective potential field provided by the lattice substrate and via which kinks get pinned due 
to lattice discreteness. Several techniques are known for this problem. According to the recent (2(j)(the collective- 
coordinate method), lattice discreteness act on the kink in two different aspects namely on the kink trajectory (time 
dependence of the kink center-of-mass coordinate within the discrete PNP), and shape. The renormalization of the 
kink shape(discreteness effect-induced kink shape dressing) is though to be manifest of radiations by the pinned kink, 
of discrete phonons as a shape response to the lattice discrete structure. To our knowledge there is yet a sound 
experimental evidence of such phenomenon, so we will not consider this aspect. Moreover, we are interested in the 
kink static configuration which is a suitable approximation of the " pinned-kink" configuration so regular in structural 
transition processes. Thus we are left with the sole problem of calculating the PNP. We start by writing the total 
energy of the discrete system considering the static part of eq. (Q): 



C7 = ^G(n), G{n)=M Cl(^-\ 
„— 1 \ / 



The discrete sum in (|3|) can be evaluated using the Poisson sommation formula J_2|,|_3| , 

U = M C 2 jr j dn (^r 1 ) exp (2mmn) 

m — — 00 00 ^ ' 

The integrand in this last equation is derived from the single-kink solution eq. (^|), which gives: 
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#»V \ B o 1 R 2 21(1+ /J 2 ) , m , 

B o = 1T7-T— ( 10 ) 



dn ) \ cosh(^) + cosh(X ) J " d2 iv) 

We take only the real part of eq. (||) which is equivalent to restrict the sum over cosine harmonics. Otherwise, this 
is consistent with the fact that G(n) is an even function of n. It turns out that the coefficients of this series must be 
determined by the following class of integrals: 
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co S (XY m ) [7 sinh 3 / 2 pf Q )r(2 + zY m )T(2 - zY m ) _ 3/2 
dX [ COS HX) + coMX )? = V 2 f(2) P_ 1/2+jy Jcoth(X o) ] (11) 
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where p is the Associate Legendre function |2l| . The complex-argument Gamma function Y appearing in this set of 
integrals is defined such that: 

r(2 + iY m )T(2 - iY m ) = 1 r F(2) |2 y2 1 (13) 
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The product function is evaluated by the standard formula p!| ; 
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Rewriting (^) as a Fourier series in m, combinations of eqs. (Jl0|)-(|l4j) lead to the following expression of the Fourier 
components: 



2d(/i) 0? sinh 3 / 2 (X )B 2 M C 2 



u ™ = ^ \ ^^Sf^ V 5*=%^ [coth(X ° )] (15) 

The most indicated framework where to check the consistency of our analytical theory is the numerical treatment. 
We performed the integration in (^) following an extended midpoint scheme. Results are displayed on figures 1 and 
2 for arbitrary values of the deformability parameter \i. The main variable here is the dimensionless kink bare width 
=¥■. dk corresponds to the asymptotic limit (i — > 0, where the two models reduce to the well-known cf> 4 model. Figures 
la and lb are the kink continuum energies at rest i.e. = Uo, while figures 2a and 2b represent the PNP barriers 
Epn = U\. Curves are plotted in unit of M C C 2 . Though the figures indicate rather complex behaviours of the PNP 
barrier with respect to /i, their most stricking feature is the drastic exponential fall-off of the pinning barrier with 
increasing fj,. For the first model, this feature relates to the fact that the greater the deformability parameter the 
narrower the kink. On the contrary, the kink width in the second model is independent on /i but however, one clearly 
sees that the corresponding PNP barriers show almost the same dependence on fi as in the first model. 
The drastic fall-off of the PNP barriers can be interpreted in terms of the kinks the overcoming lattice pinning effect. 
Such a behaviour suggests an enhancement of the conducting regime in the discrete system at relatively narrow kink 
widths. This sounds consistent with previous predictions |T^,|l^,|l3|,^0| within the 4> 4 theory, otherwise, by tuning 
the deformability parameter (which in the context of ID molecular systems may correspond to adjust the equilibrium 
positions of ions, for instance), it is expected that the correct order of magnitude of transport parameters of DWP 
systems will be obtained with relatively best accuracy p3|. 
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I. FIGURE CAPTIONS: 

Figure la: Reduced continuum energy of a kink at rest for the first PDWP model in the text. 
Figure lb: Reduced continuum energy of a kink at rest for the second PDWP model. 
Figure 2a: The reduced Peierls-Nabarro potential amplitudes for the first PDWP model, 
figure 2a: The reduced Peierls-Nabarro potential amplitudes for the second PDWP model. 
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